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I. INTRODUCTION 


The discovery of the cosmic acceleration 0 has gener¬ 
ated a lot of interest in Friedmann cosmological models 
with various kinds of exotic forms of matter, collectively 
known as dark energy [2]. Unsurprisingly, some of these 
models allow for novel type of singularities (occurring in 
finite time), which can be either strong or weak 1, 0, 
1 , 0 . 

In particular, the Big Rip singularity occurs at infinite 
scale factor, diverging Hubble parameter H and diverg¬ 
ing H 0, i- This singularity is strong according to 
both Tipler’s 0 and Krolak’s 0 definition. Some other 
singularities are alike, but they occur at finite scale fac¬ 
tor (FSF). The FSF singularity 0, 0 is also strong ac¬ 
cording to Krolak’s definition, but is classified as weak 
according to Tipler’s definition. 

The sudden future singularity (SFS), discovered many 
years ago 0, is characterized by both a finite scale 
factor and a finite Hubble parameter, while only H di¬ 
verges 0 - As consequence the energy density is finite, 
while the pressure diverges at SFS. The geodesic equation 
(containing only H) remains regular at SFS, such that 
point particles will pass through the singularity. How¬ 
ever for an infinitesimally short time the tidal forces in 
the geodesic deviation equation become infinite due to 
the diverging H. The SFS singularity is weak in both 
Tipler’s and Krolak’s definitions, hence it has been con¬ 
jectured that finite object may also cross this singularity 

A subclass of the SFS singularities is the Big Brake, 
introduced in Ref. fl2| as emerging in the presence of 
a particular tachyonic scalar field (i.e. Born-Infeld-type 
field). Such a singularity is characterized by a full stop 
in the expansion, occurring at finite scale factor and van¬ 
ishing energy density, augmented by a diverging deceler¬ 
ation and pressure. There is also the w-singularity Q, 
0 , when both the pressure and energy density vanish, 


however their ratio, the barotropic index w diverges, is 
also weak in both definitions. 

Interesting features of the SFS singularities have been 
found on the example of the flat Friedmann universe filled 
with a combination of cold dark matter (dust) and tachy¬ 
onic field, the latter obeying a particular potential of 
trigonometric functions 1121 ]. The model possesses two 
surprising features. First, for certain model parameter 
range, at some stage of the cosmological evolution, the 
tachyonic field is constrained to transform into another 
Born-Infeld-type field, the pseudo-tachyon. Secondly, in 
this regime the universe unavoidably evolves into the 
SFS. 


In this model in the absence of the dust component 
the SFS singularity represents a Big Brake. Confronta¬ 
tion with type la Supernovae (SNIa) data showed that 
the allowed cosmological evolutions can reach either the 
Big Brake or the de Sitter attractor 0. Further, in 
Ref. [10 it was proven that the Big Brake singularity is 
traversable, and the evolution continues through a col¬ 
lapsing universe into a Big Crunch. 

When the dust component is included, the SFS can¬ 
not be a Big Brake anymore 0. Indeed, the universe 
reaches the singularity with a finite but nonzero H. Its 
traversability is guaranteed by certain changes in the 
matter properties [13] , [10 , imposed by the requirement 
of a smooth passage through the SFS. These changes 
lead to the notion of the quasi-tachyonic Born-Infeld- 
type field, a phantom field [10 with negative energy den¬ 
sity. A throughout confrontation of these models with 
SNIa, baryonic acoustic peaks, and cosmic microwave 
background temperature anisotropy power spectrum did 
not disrule the evolutions leading to either the de Sitter 
attractor or to the SFS [10 . 

The present work represents a continuation of the se¬ 
ries of papers 0, 0-0, 0, investigating the conse¬ 
quences of the introduction of a tiny positive curvature, 
which remains compatible with the recent Planck data 
21j, 1221. As will be shown, the addition of the curva- 
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ture term enriches the velocity phase space of the possi¬ 
ble cosmological evolutions, adding to it regimes which 
were impossible to realize in the flat model. We study 
these regimes numerically and show the existence of a 
critical behavior in the formation of the singularities in 
the framework of these tachyonic cosmological models. 
Such critical phenomena were well-known in the numer¬ 
ical study of gravitational collapse of a scalar field [23] 
and of various other fields (as reviewed in Ref. [2J|), 
with the critical evolutions separating black-hole reach¬ 
ing evolutions from those which do not lead to black hole 
formation. 

The structure of the paper is as follows. In Section 
Ull we introduce the equations of the tachyonic cosmol¬ 
ogy in the presence of dust and curvature. In Section 
m we confront the luminosity distance-redshift relation 
with the supernova data. In Section m we investigate 
numerically the dynamics and give the classification of 
all possible regimes, with emphasis on the new types of 
evolutions. The results are summarized in the Conclu¬ 
sions. 

The units are chosen as c = 1 and 8irG/3 = 1. 


II. A MIXTURE OF DUST AND TACHYONIC 
FIELD IN A FRIEDMANN UNIVERSE 

We consider a Friedmann universe 

ds 2 = dt 2 —a 2 ( t ) — ^J^dr 2 + r2 d0 2 + r 2 sin 2 ddq , 

(!) 

with cosmological time t, comoving radial distance r, po¬ 
lar and azimuthal angles 6 and </>, unnormalized curva¬ 
ture index I\ > 0,= 0 or < 0 and scale factor a(t). The 
combined energy momentum tensor is an ideal fluid with 
energy density p and pressure p. The dynamics of the 
scale factor is governed by the Raychaudhuri equation: 

• K 3 , 

H-—^ + -(p + p)—0, (2) 

a z 

and the change in the energy density by the continuity 
equation: 

p + 3H(p+p) =0 . (3) 

The overdot denotes time derivative and H = a/a is the 
Hubble parameter. The first integral of the system ([2]) - 
© is the Friedmann equation: 

H 2 -\ —- = p . (4) 

a z 

We assume that the Friedmann universe is filled with a 
mixture of tachyonic field and cold dark matter (CDM). 
The mixture of the two ideal fluid components (dis¬ 
tinguished by m and T subscripts) is characterized by 
P = Pm + Pt and p = pt, where 

< 5 > 


with integration constant p* m . The evolution of the homo¬ 
geneous tachyonic field T(t ) follows from the Lagrangian 
density 125!]] 


L T = -V=gV(T) sjx-g^s 2 , (6) 

with s = T. Variation with respect to the metric gives 
its energy density 


Pt 


ViT) 


( 7 ) 


and pressure 


p T = -V (T) Vl-s 2 . (8) 

As w = p/p = s 2 — 1, in principle the tachyonic field can 
mimic dark energy. The particular case of a trigonomet¬ 
ric tachyonic potential [12 ] 


with 


Ay/1 - (1 + k)y 2 

1 -y 2 


y = cos 


^A(l + k)T 


( 9 ) 


( 10 ) 


A > 0 and — 1 < k < 1 model parameters further leads 
to interesting behavior, including a smooth passage into 
the s > 1 regime and a sudden future singularity which 
can be passed through by geodesics which then bounce 
back into a recollapse [12]-j20|. 

Note that the system is invariant under the simulta¬ 
neous changes (y —> — y, s —> — s).This yields a double 
coverage of the dynamics in these velocity phase space 
variables. 


III. TEST WITH SUPERNOVAE DATA 


The evolution of the Friedmann universe is determined 
by the initial conditions H 0 = H (z = 0), yo = y (z = 0), 
so = s (z = 0), (z is the redshift), the present values of 
the cosmological parameters 


Htv. — 


aim 


Q k = — 


K 

aim 


( 11 ) 


associated to the CDM and the curvature, and the model 
parameters A, k. The parameter Qk is constrained by 
observations to 95% confidence level: —0.0065 < fix < 
0.0012 by WMAP data [13 and -0.019 < fl K < 0.011 
by Planck data (see fifth column in Table 5 of [22|). In 
what follows we set f l m = 0.315 (supported by the cosmic 
microwave anisotropies measured by the Planck satellite 
[Uli H3)i ^ = 0-44, H 0 = 70 km/sec/Mpc and we ana¬ 
lyze both Qk = —0.0065 and 0 (the former being com¬ 
patible with K > 0, the latter representing a spatially flat 
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universe). The Friedmann equation ([4]) represents a con¬ 
straint among y 0 , s 0 and A, thus only the first two are 
varied. The confrontation of the tachyonic model with 
the Union 2.1 SNIa data set [26| through a y 2 -test as 
described in Ref. [l5| relies on fitting with the computed 
luminosity distance (dL)-redshift relation. A dimension¬ 
less luminosity distance dp = H^l is given by 


= x(z) for K = 0, (12) 

1 + z 

\/-n K = sin (x (z)) for K > 0 , (13) 

1 + z 

with 

dx (z) = V-^k , 14 x 

dz H 

(H = H/Hq) We show the fit on Fig |T] The contours 
refer to the 68.3% (ler) and 95.4% (2cr) confidence levels. 



FIG. 1: (Color online) The fit of the luminosity distance vs. 
redshift for the parameters k = 0.44 and fi m = 0.315 in the 
parameter plane ( yo, w 0 = (l + so) ) (upper panel for 11 k = 
—0.0065, lower panel for VLk = 0). The contours refer to the 
la and 2 <t confidence levels. The colour code of y 2 is indicated 
on the vertical stripes. 


Thus the introduction of the curvature does not ham¬ 
per the compatibility of the model with the SNIa obser¬ 
vations. 


IV. VELOCITY PHASE SPACE 


The y 2 < 1/(1 + k) , s 2 < 1 domain of the velocity 
phase space represents proper tachyonic fields, character¬ 
ized by the Lagrangian ([6]), energy density © , pressure 
(|g|) and potential ©• A subset of the initial conditions 
yo, so allow the universe to evolve outside this regime 
through the corner points y, s = ±1 where the potential 
V vanishes, as described in Ref. [jjj]. These are excep¬ 
tional points, where the curvature does not diverge when 
s 2 = 1. When s 2 > 1, (1 + fc)” 1 < y 2 < 1 the field 
becomes pseudo-tachyonic and V imaginary. Although 
the Lagrangian © does not change in this regime, it is 
convenient to introduce the real potential 

AW(1 + k) y 2 — 1 

W = iV=^l -- . i5 

i - r 

Note that in the Lagrangian W absorbs the imaginary 
factor arising from the change of sign in the equation of 
state w = s 2 — 1. The expressions of the energy density 
and pressure are also unchanged in the pseudo-tachyonic 
regime. Rewriting the Lagrangian, the energy density 
and pressure in terms of W however manifestly shows 
that all these quantities stay real: 


Lp = yjg u s 2 - 1 , 

W(T) 

pp = 7 ^= 1 ’ 

p P = W (T) y/s 2 - 1 . 


(16) 

(17) 

(18) 


The pressure is positive, thus the expansion of the uni¬ 
verse slows down in this regime. 

Further evolution to s —> ±oo leads to SFS, arising for 
finite values of 2 , H and W and pp —> 0, pp —> oo. As 
was shown in Ref. [18] the evolution can be continued 
smoothly across these soft singularities at the price of a 
sudden change of the Lagrangian, which will represent a 
quasi-tachyonic field: 


l Q = \[ ZI gW\Jg tt s 2 + 1 , 


PQ 


-W{T) 
Vs 2 + 1 ’ 


PQ = W (T) Vs 2 +1 . 


(19) 

( 20 ) 
( 21 ) 


Note that wq = Pq/pq = — (s 2 + l) / w, thus the 
quasi-tachyon is a phantom with negative energy density. 

In order to formulate the dynamical equations with 
dimensionless variables, we normalize the potentials as 


v, _ - (1 + k)y 2 

1 -y 2 
+ k)y 2 - 1 
1 -y 2 

(i.e. n A = A /H 2 , V = V/H 2 , w = W/H 2 ). 


( 22 ) 

( 23 ) 
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The dimensionless variables z(t) and H(t) evolve as 

(] Z /s. 

- s + (l + z)H = 0, (24) 

at 

“T' - + 9 {pm + P +p) + (1 + z) = 0, (25) 

at 2 

with t = Hot. Eq. (HH1) follows from the definition of 
z, while Eq. (1251) is the Raychaudhuri equation. Here 
p m = I2 m (1 + z) is the normalized density of the dust, 
while p stands for either of the normalized densities 


Pt — 

V 

(26) 

71 - s 2 ’ 

Pp = 

w 

(27) 

a/s 2 - 1 

PQ = 

—W 

(28) 

a/s 2 + 1 


holding for the tachyonic, pseudo-tachyonic or quasi- 
tachyonic regimes, respectively. Similarly, p denotes one 
of the normalized pressures 


Pt 

= -77i-s 2 , 

(29) 

Pp 

= l77s 2 - 1 , 

(30) 

PQ 

= l77s 2 + 1 . 

(31) 


The dynamics of the tachyonic field y(t) emerges from 
the definition of s as 



(1) (2) (3) (4) (5) (6) (7) 



y 


+ ^s\/^a (1 + k) (1 — y 2 ) — 0 , (32) 

while the evolution of s(t) in the tachyonic, pseudo- 
tachyonic and quasi-tachyonic regimes are given by the 
following Euler-Lagrange equations: 

| + (l- s2 )( 3 ^ + f) = °> ^ 

S +(i_s2) ( 3ffs+ ^) = 0; (34) 

S +(i+s2) ( 3ffs '^) = °- (35) 

In the neighborhood of the soft singularity y can be ap¬ 
proximated as 


yp = 


VQ = 


, 76 r -j 

ys + — v 1 - ys 


x v^A (1 + k)j-—Vts - t , (36) 

V Hs 


ys 


76 


- 7T'/ l ~ys 2 


X-y/0 A (1 + k)\ — y/t - t S . 

J~LQ 


( 37 ) 


FIG. 2: (Color online) Upper panel: The numerical evolution 
of the dynamical equations yields seven type of trajectories 
(I, Ila, lib, Ilia, Illb, IV and V), separated by separatrices 
(red). Lower panel: the pseudo-tachyonic evolutions of the 
upper right corner. We chose Oa = 0.8. 


(These formulas are analogous to Eqs. (30) and (50) of 
Ref. [IH presented in terms of the original tachyonic 
variable T, and they hold for the upper right strip of 
Fig- m where yp > ys > VQ-) The two expressions coin¬ 
cide at the singularity. In the quasi-tachyonic regime the 
Hubble variable decreases, reaching zero at some point 
which is followed by a recollapse through the same type 
of singularity to the pseudo-tachyonic regime, ending in 
a Big Crunch. This was always the case for I\ = 0 jTg] ■ 

Due to the extra gravitational attraction represented 
by 12 k < 0 however some of the pseudo-tachyonic evolu¬ 
tions do not reach the singularity, avoiding the quasi- 
tachyonic regime. Instead the recollapse into a Big 
Crunch commences earlier. These two types of evolu¬ 
tions are separated by a critical trajectory when the soft 
singularity is reached exactly with H = 0. In what fol¬ 
lows we will investigate this criticality. 

The numerical evolution of the dynamical equations 
yields seven types of trajectories, I, Ila, lib, Ilia, Illb, 
IV and V, represented on the upper panel of Fig. [2] 
The lower panel magnifies out the pseudo-tachyonic evo- 
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FIG. 3: (Color online) The trajectories presented on the 
lower panel of Fig. [2] are extended to infinity, and mapped 
into the interval [0,7r/2] by a conformal transformation a = 
arctan(10~ 3 s). The pseudo-tachyonic and quasi-tachyonic 
regimes of each of type IV, Ilia, Illb and V evolutions (blue), 
and the separatrices between types IV-IIIa, Illa-IIIb and 
Illb-V (red) are shown. The soft singularities correspond to 
a = 7 t / 2 . The evolutions labeled by numbers are identical to 
those represented on the lower panel of Fig. [2] The critical 
evolutions leading to Big Brake singularities in the presence 
of dust are (4) and (6). 

lutions of the upper right corner. 

The traject ories of t ype I originate from the Big Bangs 
{y,s) = (1, yj 1 + 1 /k) and (y,s) = (-1 ,-y/l + l/k) in 
pseudo-tachyonic regime, then evolve into the tachyonic 
regime towards a de Sitter expansion, represented by the 
origin (y,s) = (0,0). The trajectories of type Ha orig¬ 
inate from the lines s = ±1 representing Big Bang sin¬ 
gularities, then evolve from the pseudo-tachyonic regime, 
through the corner points into the tachyonic era, where 
they also reach the de Sitter attractor. The trajecto¬ 
ries of type lib are similar, but they never cross into the 
pseudo-tachyonic regime 

The trajectories of type Ilia and Illb also originate 
in Big Bangs at the lines s = ±1, pass through the 
tachyonic regime and through the cornerpoint into a 
pseudo-tachyonic regime, followed by a collapse into a 
Big Crunch (on the y = ±1 line). The difference be¬ 
tween them is shown on Fig. [3] While Ilia reaches the 
soft singularity, evolves through a short quasi-tachyonic 
regime in which the universe start to recollapse towards 
a second soft singularity, those of type Illb never reach 
a soft singularity. 

The trajectories of type IV originate in a Big Bang 
located in the pseudo-tachyonic regime (y = ±1 line), 
evolve into the tachyonic regime, then return into a 
pseudo-tachyonic regime, pass twice through soft singu¬ 
larities (with a short quasi- tachyonic period), finally col¬ 
lapse into a Big Crunch. The trajectories of type V are 


born in the Big Bang singularities on the lines s = ±1, 
after which they mimic the evolution of the trajectories 
of type IV following their second passage into a pseudo- 
tachyonic regime. 

On Fig. [3] we represent all types of trajectories which 
will reach a pseudo-tachyonic regime in the future (types 
Ilia, Illb, IV and V). From among these Ilia, IV, V 
pass twice through soft singularities, with an intercalated 
quasi-tachyonic regime, while those of type Illb never 
reach a soft singularity. We also determined numerically 
the critical evolutions separating these two types of be¬ 
haviors. They are the evolutions leading to Big Brake 
type singularities in the presence of dust. 


V. CONCLUDING REMARKS 

In this paper we have expanded previous studies of a 
tachyonic cosmological model with trigonometric poten¬ 
tial and dust exhibiting a rather rich spectrum of possible 
evolutions, by including curvature. As expected, a posi¬ 
tive curvature could counterbalance the influence of the 
dust (manifesting itself in a positive H and nonvanishing 
energy density at the SFS), in the sense that it can com¬ 
pensate for the nonzero energy density at the SFS. When 
the curvature term is strong enough, it stops the expan¬ 
sion before the encounter with the SFS, preventing its 
formation. If the curvature term is less strong, then the 
cosmological SFS is still formed, followed by the transi¬ 
tion to the quasi-tachyonic regime, but its slowing-down 
effect is expected to reduce the duration of this regime 
as compared to the flat case. 

We have shown by numerical methods that it is possi¬ 
ble to have a particular balance between curvature and 
dust effects, leading to critical evolutions, represented by 
the curves (4) and (6) on Fig. [3] When such a criti¬ 
cal evolution is realized, the universe reaches exactly a 
Big Brake singularity, which was forbidden in the flat 
case. The critical evolutions reaching a traversable Big 
Brake separate evolutions into SFS from completely reg¬ 
ular ones. In this sense they resemble the criticality of 
singularity-forming in gravitational collapse. This criti¬ 
cality is a new feature of tachyonic cosmologies, induced 
by curvature. The possibility to reach a Big Brake type 
singularity when dust is present is the second new fea¬ 
ture of the closed Friedmann universes with the tachyonic 
field. 
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